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Abstract 

We study the likelihood of energy condition violations in the history of the Universe. Our method is 
based on a set of functions that characterize energy condition violation. FLRW cosmological models 
are built around these "indication functions" . By computing the Fisher matrix of model parameters 
using type la supernova and Hubble parameter data, we extract the principal modes of these func- 
tions' rcdshift evolution. These modes allow us to obtain general reconstructions of energy condition 
violation history independent of the dark energy model. We find that the data suggest a history 
of strong energy condition violation, but the null and dominant energy conditions are likely to be 
fulfilled. Implications for dark energy models are discussed. 

Keywords: cosmology: theory - dark energy - gravitation - rclativistic processes 



1. INTRODUCTION 

The energy conditions introduced by lHawking fc EUii 
(|1973l Chapter 4) are coordinate-invariant inequality 
constraints on the energy-momentum tensor that ap- 
pears in the Einstein field equation. They play an essen- 
tial role in the proof of theorems in the general theory 
of relativity c oncerning th e existence of spacetime sin- 
gularities (see iWakl 119841 ). Due to their simplicity and 
model-independence, the energy conditions are listed as 
one of many approaches to understand the evolution of 
universe. 

Of the many proposed energy conditions, the ones we 
analyze in this paper are the null, the strong, and the 
dominant energy conditions (abbreviated respectively as 
NEC, SEC, and DE C). These ene rgy conditions can be 
expressed as follows (|Carrollll2004 Chapter 4) for a gen- 
eral stress-energy tensor T: 

• NEC: T a/3 n a n f:S > for all null vectors n; 

• DEC: T af3 t a t > and T a aTh a V < for all time- 



i/3 

like vectors t; 



SEC: Tap^t 13 > \T^t s t s for all timelike vectors t. 



The energy conditions are of great relevancy to the study 
of cosmology. For example, the NEC is an important 
condition of stabi lity for fluids (|Dubovskv et al.l 120061 : 
iBuniv et al.l I2006D . hence it is useful for the analysis 
of cosmological models. The DEC guarantees stability 
of a source obeying it and imposes on the dark energy 
equation of state p arameter w a lower bound w > — 1 
([Carroll et al.ll2003| ). Violation of the w > -1 bound 
can le ad to the "big rip do omsday" scenario of the Uni- 
verse (ICaldwell et all 12001 . It can also be shown that 
sudden future singularity solutions of the Friedmann- 
Lcinaitrc-Robcrtson- Walker (FLRW) cosmology always 
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violate the DEC (|Lakell2004| ). Finally, SEC violation is a 
typical tra it of a positive cosmological constant A (for ex- 
amp le. seelLi et al.l 1201 If) and other dark energy models 
(see iSchuecker et al.l 120031 and references therein) . Wc 
therefore find an investigation of these energy conditions 
useful for our understanding of the Universe's evolution. 

The knowledge of cosmic energy condition violation 
has be e n great ly facilitated by precise observational data. 
Visscr (1997a. b) shows that Hubble constant and stellar 
age measurements, which bound the age of the Universe, 
suggest a history of SEC v iolation (see also the review by 
IKrauss fc Chab ovcr 20Q3). X -ray galaxy cluster number 
count and type la supernova (SNIa) luminosity distance 
data have been applied to a constant-w model, which 
also suggest SEC vi olation without signific ant indication 
for NEC violation (ISchuecker et all 120031) Similar re- 
sults are also obtained bv lLima et al.l (I2008al) using avail- 
able SNIa data at that time. On the other hand, if one 
assumes the energy conditions, constraints on a variety 
of cosmological observablcs or parameters can be pre- 
dicted, such as the Hubble parameter, luminosity and 
angular diameter distances, lookback time, total density 
parameter f2(z), energy density p(z), and pressure p(z) 
(jCattoen fc Visserl[200a i. 

In this paper we study the issue of energy condition vi- 
olations from a data-driven perspective. Specifically, we 
aim at reconstructing what we can say about the likeli- 
hood of energy condition fulfillment (or violation) given 
the data. To achieve this, we use the energy conditions 
themselves to construct a family of descriptive cosmologi- 
cal models for the recent history of cosmic expansion, and 
subject them to statistical test. Throughout our paper 
we assume a universe of perfect fluids that lead to the 
spacetime solution characterized by the FLRW metric. 
We use t he Union2 SNIa lumino sity distance data com- 
pilation (|Amanullah et all 120101 ) and the observational 
Hubble parameter data from different i al ages of passively 
evolving galaxies (jSimon et al.1 120051 iStern et "aD I2010T ) 
and ra dial baryon acoustic o scillation (BAO) measure- 
ments (|Gaztanaga et al1l2009l ). The cosmological param- 
eters used in this paper: f2 m = 0.274, = 0.726, and 
Hubble constant H = 74.2 km s" 1 Mpc -1 

This paper is organized as follows. In Section [2] wc 
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specify our models built from our indication functions of 
energy condition violation and provide solutions for the 
cosmic expansion rate. We proceed to Section [3] where 
we lay out the procedures of analysis using luminosity 
distance and Hubble parameter data. Our main results 
are presented in Section 2] and are discussed in Section 

El 

2. ENERGY CONDITIONS AND COSMOLOGICAL 
MODELS 

2.1. Overview 

As summarized by ICarrolU ((20041 Chapter 4) , the en- 
ergy conditions can be expressed by simple inequalities 
in terms of the energy density p and the pressure p in the 
setting of a homogeneous and isotropic FLRW universe: 

• NEC: p + p > 0, 

• DEC: p > and -p<p< p, 

• SEC: p + 3p > and p + p > 0. 

By virtue of the Fricdmann equation, the energy con- 
dition bounds can be expressed in terms of the Hub- 
ble constant-normalized, dimensionless expansion rate 
E(z) — H(z)/Hq and its first derivative, or the decel- 
eration parameter defined by 



• NEC: 



1 + z dE 
E\z) dz 



(1) 



The results, following ILima et al.l (|2008aD . are given be- 
low: 



NEC: 



DEC: 



, , n k (i + z) 2 , „ 

9 <*>-^W- +1 - 0> 



2Q k (l + *) 2 . 
2 - c j( z ) T^J-^ > °> 



(2) 



(3) 



, , n k (i + z) 2 

F NE c{z) = q(z) h 1 ' 



(5) 



DEC: 



F DEC (z)=2-q(z) k V,. , ^ , (6) 



E 2 (z) ' 



• SEC: 



Fsec(z) = q{z). 



(7) 



With the introduction of these indication functions, 
inequalities ([2](4j) become ordinary differential equations 
that can be integrated to find E(z). We may look at 
them as mappings from the set of dimensionless Hub- 
ble expansion rates E(z), to that of indication functions 
F(z), and vice versa. These mappings enable us to re- 
construct F(z) from observables related to E(z). 

Our formalism so far remains general enough for any 
model that makes a prediction on E(z). In particular, 
we will concentrate on the underlying evolution of E(z) 
that is favored by observational data0- This approach 
has bee n used on the study of possib l e evolution of dark 
energy (|Huterer fe Starkmanl 120031: IHuterer fc Cooravl 
120051 ) and the reconstruction of histori cal deceleration 
parameter q(z) ([Shapiro fc Turnerll2006l ). 

Following this line of inquiry, we can use a similar 
method to reconstruct the indication functions F(z), 
where the function to be reconstructed is "coarse- 
grained" as piecewise-constant. While this can be done 
using more elaborate function s such as the smooth , 
hyperbolic tangent function ([Crittenden et al.1 120091 : 
IZhao fc Zhanll l2010h . we find the piecewise-constant 
functions greatly simplify the integration of equations ((5} 
[7]). In the interval in which F(z) is constant, F(z) = r, 
the general solutions are given by 

• NEC: 



C{1 + z) 2r_ n ^ + *)\ (8) 



• SEC: 

q(z) > 0. (4) 

A notable feature of these bounds is that they are only 
explicitly dependent on one arbitrary constant, namely 
the curvature parameter Slk, although the matter content 
of the FLRW universe does control the functional form of 
E(z). However, for our purpose we do not need to make 
assumptions on the constituent matter of the universe 
being modeled. 

2.2. Models and Their Solutions for E{z) 

We are therefore motivated to introduce the "indica- 
tion functions" F(z) that quantifies whether an energy 
condition has been violated. Fulfillment (or violation) 
of an energy condition should be indicated by F(z) > 
(or F(z) < 0). The simplest choices are just the left- 
hand sides of inequalities ([2]|4]). We denote them as 
-^NEC) -FdeC; and -Fsec respectively: 



• DEC: 

£decW - ] /c(l + z)m-r) + m ^ 1 _ + r Z)2 , (9) 

• SEC: 

E SEC (z) = ^C(l + z)W) 7 (10) 

where C stands for the integration constant. To obtain 
the special solutions of E(z) for the piecewise-constant 
F(z) models extending from z = to an arbitrary red- 
shift, we can recursively apply the above solutions and 
the continuity condition of E{z) across subinterval end- 
points starting with E(z = 0) = 1, thereby fixing the 
integration constants. 

4 The practice of deriving a model of cosmic expansion by phe- 
nomcnologically treating the model parameters without physical 
prescriptions is widely known as the "model-independent" ap- 
proach, which does not actually exclude the use of a model. 
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3. ANALYSIS 

As we have briefly stated in the preceding section, we 
use the Fisher matrix of the full array of model pa- 
rameters to extract a small number of principal com- 
ponents that preserve the information contained in the 
data without introducing serious over-parameterization. 
The Fisher matrix elements for the model parameters 
(not to be confused with the symbol for the indication 
functions F) arc expressed by 



F=(-^) 



(ii) 



where P is the posterior probability density function 
(PDF) in the parameter space and Oi is the ith model 
parameter, and the angled brackets stand for statistical 
averaging. The principal components corresponds to the 
eigenvectors of the Fisher matrix. The posterior proba- 
bility P is to be found from observational data. 

3.1. Posterior Probability 

To construct the posterior probability P from obser- 
vational data, we assume that the uncertainties assigned 
to the measurement results are Gaussian. Under this as- 
sumption, the posterior probability can be expressed by 
an additive x 2 statistic with x 2 — —2 In P. The generic 
form of x 2 under this assumption can be laid out as fol- 
lows: 



X 



E 



th 



AT 



bs\2 



(12) 



where the symbols X th,obs denotes, respectively, the the- 
oretical or observed value of the ith observable (direct or 
indirect), and Oi is the ith Gaussian variance or uncer- 
tainty given by the data. The summation is done over 
all individual data entries. 

In dealing with parameterized models, it is natural 
to introduce nuisance parameters to be eliminated later 
by marginalization. The nuisance parameters and their 
marginalization distorts the simple form of equation (|12p 
as well as introducing computational complexities that 
often calls for Monte-Carlo techniques. Nevertheless, be- 
fore engaging any numerical computations we can elimi- 
nate certain nuisance parameters by analytical marginal- 
ization, as illustrated in the following sections. 

3.1.1. Luminosity Distance 

The luminosity distance data, such as the Union2 
dataset that we use in this paper, are usually presented 
as tabulated distance moduli with errors. Physically, the 
luminosity distance modulus is the difference of apparent 
magnitude m and the absolute magnitude M, 



fj,(z) = m — M = — 5 + 51gdi(z), 



(13) 



with oIl measured in the units of 10 parsecs. It can be 
related to the modeled expansion rate of the universe, 
E(z), by the formula 



d L (z) 



1 + z 



H 



dz' 



E(z') 



(14) 



where the sinn function is a shorthand for the definition 

'sinh.T, Ok > 0, 
sinn(a;) = i x, Ok = 0, (15) 



sin a;, 



o k = o, 

O k < 0. 



For notational convenience we employ a variable 



m(z) = 5 lg 



l + z 



: sum 



dz' 



(16) 



where the parametric dependency on the Hubble con- 
stant Hq has been separated out, leaving a dimcnsionlcss 
quantity as our model prediction. With this definition 
the modeled distance modulus can be expressed as 



/i th (z) = m(z) + 5lg H + M . 



(17) 



In this expression we introduce a numerical constant Mq 
to accommodate the numerical constants arising from 
unit conversions, as well as an additive variable char- 
acterizing the uncertain variability or spread of the stan- 
dard candles' absolute magnitudes. By doing so, we in- 
clude Mo as another nuisance parameter which affects 
the overall scaling of the modeled cosmological distance. 

The x 2 statistic for a luminosity distance modulus 
dataset, according to equation (fT2|) . is therefore 



X 



E 



[m{zi) + 51gH + M 



obs 



(*)] ! 



(18) 



The above expression is quadratic in the combination of 
nuisance parameters (51g\ffo + Mo): 

X 2 = A(5 lg H + Mo) 2 + 25(5 lg H + M ) + C, (19) 

where, for notational simplicity, we have defined the sym- 
bols 

^ = E^> 

m{zi) - ii ohs { Zl ) 



B 



E 



(20) 



c = E 



[m{zi)-ti ohs {zi)\ 



The nuisance parameters Ho and Mo can now 
be marginalized over by Gaussian integration over 
(—oo, oo), leaving a "reduced x 2 " in the form of 



x 2 



c 



A ' 



(21) 



This expression of x i s a general result for lu- 
minosity distance data (Idi Pietro fc Claeskensl 120031 : 
iNesseris fc Perivolaropoulosll2004D . 

3.1.2. Hubble Parameter 

The observational Hubble parameter data we use 
are t abulated in Refs. (|Zhang et all 120101 : IMa fc Zhangl 
20111) and form an independent dataset. As noted by 



Zhang fc Zhul (|2008l ). the Hubble parameter data could 



help us catch possible "wiggles" of E{z) predicted by cer- 
tain dark energy models that would have been flattened 
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Figure 1. Binning of the rcdshift range to accommodate enough 
data in each of the bins. The histogram shows the number of SNIa 
data points distributed across the redshift range, with each bin 
covering width 0.1. 

out by the luminosity distance test. For the Hubble pa- 
rameter data our model requires equation (|12[) be in the 
form of 



x ^ [HpEjzj) - H obB (zi)Y 
2 - 



(22) 



We adopt the result of H = 74.2 ± 3.6 km s" 1 Mpc -1 
measured from nearby Cep heids by the Hubble Space 
Telescope (jRiess et al.ll2009t) as a Gaussian prio r and use 
the analytic expression given by IMa fc Zhang] (|2011D to 
construct the \ 2 with Hq marginalized over. This statis- 
tic from H(z) data is added to the one from SNIa data 
(eq. [21]) to give the final posterior log-PDF. 

3.2. Redshift Binning 

Ideally, the coarse-graining approximation of the F(z) 
evolution should approach continuity as its limit. How- 
ever, for our analysis using the data of 557 SNela from 
Union2 and the 13 H(z) measurements, the redshift dis- 
tribution of data poses a limit on how fine we can dissect 
the redshift range into bins, hence the resolution of our 
reconstruction of F(z). Each bin should cover enough 
data to make the value of F(z) well-constrained inside 
of it. To this end, we introduce 11 bins equally spaced 
to cover the redshift range, shown in Figure [T] The first 
10 bins are kept to be equal in width, but the last one is 
widened to cover the diminishing tail of the SNIa redshift 
distribution. 

3.3. Evaluation of the Fisher Matrix 

The Fisher matrix elements as defined by equation (|11[1 
are statistical averages involving unknown parameters. 
However, for our purpose we cannot evaluate them by 
averaging in the absence of a known posterior. In fact, 
it is neither practical nor logical to do this, as our final 
models to assess the evaluation of energy conditions are 
not the model described so far per se, but the models to 
be built from the principal components extracted from 
the Fisher matrix. 

In practice, the Fisher matrix elements can be esti- 
mated at point in the parameter space where the pos- 
terior PDF is close to the maximum. In other words, 



we can approximate the numerical values of the Fisher 
matrix elements using a fiducial model, which we choose 
to be the concordance flat ACDM model with Q m = 
0.274, Oa = 0.726 as favored by combined observat ional 
data (jAmanullah et all l20l0t iKomatsu et all 120 111) and 
H Q = 74.2 km s" 1 Mpc" 1 from (iRiess et all 120091) . In 
a similar setting, iShapiro fc Tur ner (200(|) used a even 
simpler, constant-input model as the fiducial one in their 
analysis, and found the method to be robust against the 
introducing of a simplified fiducial. Our own findings are 
in agreement to their claim. 

To add robustness and safeguard against limited sam- 
pling in the parameter space, we resampling (|Efronlll982| ) 
to check for any possible bias resulting from not treat- 
ing the averaging in equation (|llj) rigorously. The result 
confirms the robustness of Fisher matrix estimation, and 
the average from the resamplcd Fisher matrices are used 
as the "reference stack" in subsequent steps. 

The Fisher matrix elements are thus calculated using 
the x 2 (or posterior log-PDF) expressions given in the 
preceding Sections 13.1.11 and 13.1.21 Notice that the nor- 
malization constants discarded in the previous steps do 
not contribute to the matrix elements because they are 
additive constants in the log-PDF independent of the pa- 
rameters, and their derivatives vanish. 

4. RESULTS 

We numerically compute the (N + 1) x (N + 1) Fisher 
matrix for F{z) in each of the N bins and Ok- At this 
stage, the Ok parameter can be effectively marginalized 
over by projecti ng the matrix onto the N x N subspace 
(|Dodelsonll2003l Chapter 11). The eigenvectors of the re- 
sulting matrix, /j, are then retrieved. In this manner, we 
can thus build a chain of models with successively more 
parameters. The eigenvectors are the principal compo- 
nents we are after, and they are shown in Figures [2] [3] 
and [4] for three energy conditions. Then the reconstruc- 
tion we need can be expressed as a linear combination of 
the principal components 



(23) 



Notably, the first three principal components appears 
stable regardless of which energy condition the recon- 
struction is for. However, the higher-order and noisier 
modes show greater change in the form with respected to 
the energy condition being reconstructed. Furthermore, 
when reconstructing F(z) by fitting the linear expres- 
sion (Eq. [23]), the curvature parameter Ok appearing in 
equations ([5]) and ([6]) has been marginalized with numer- 
ical integration. 

The reconstructed F(z) for the three energy conditions 
are shown in Figures [5] |6] and [3 We have use only the 
first three principal components in the fitting of F(z). 
The fitted values of the coefficients are listed in Table [T] 

This cutoff from the full spectrum of principal compo- 
nents introduces a well-known artifact, namely a bias at 
the higher end of the redshift ([Shapiro fc Turned 120061 ) 
that suppresses the reconstructed value and its error bars 
towards zero. Nevertheless, we note that a value of F(z) 
close to zero is consistent with our lack of knowledge 
about the violation of energy conditions by the very con- 
struction of F[z). In other words, a vanishing value of 
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0.0 0.2 0-4 0.6 0.8 1.0 1.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



Figure 2. Four well-constrained principal components of 
Fnec( 2 )- They are normalized so that J ff(z)dz = 1. fi's (i = 
1, 2, 3, 4) are the first four eigenvectors of the Fisher matrix. 
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Figure 3. Same as Figure[2]but for -FdecC- 2 )- 
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Figure 4. Same as Figure[2]but for Fsec( 2 )- 



F(z) tells us that our bet of finding energy condition ful- 
fillment is no significantly higher or lower than that of 
finding its violation. Therefore, this bias in the value of 
F(z) does not tempt us with false claims about energy 
condition violations, but merely restates our ignorance 



Table 1 

Best-fit coefficients of the first three principal 
components 



a, 


NEC 


DEC 


SEC 


a l 


-0.51 ±0.08 


2.40 ±0.10 


0.88 ±0.08 


(12 


0.12 ±0.21 


-0.31 ±0.22 


-0.18 ±0.21 


0,3 


-0.16 ±0.45 


-0.25 ±0.50 


-0.43 ± 0.40 




Figure 5. Reconstruction of Fnec(z) using the first three eigen- 
vectors. The central blue solid line in each bin is the best three- 
mode fitting result. The gray and yellow bands show 1- and 2-cr 
uncertainties. 



Figure 6. Same as Figure |51 but for Fd^q(z). 

about it at redshift ranges where there are insufficient 
data. 

5. CONCLUSION AND DISCUSSIONS 

We have explored the likelihood of energy condition 
violation during the evolution of the universe using the 
indication functions we proposed in Section [2.1[ The in- 
dication functions are reconstructed from SNIa and H(z) 
data. 

Our result for NEC does not prefer a history of vio- 
lation, as is evident from Figure [SJ A slight but rising 
trend can be seen from z ~ 0.7 up to now. This trend 
is more pronounced in our result for DEC presented in 
Figure |6l and we can be fairly certain that there has been 
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-1.0 



Figure 7. Same as Figure [5] but for Fsec( 2 )- 

no DEC violation for z < 0.5. 

However, the result for SEC shows fairly strong indica- 
tion of violation. This shows the view of an accelerating 
cosmic expansion is compatible with our combination of 
data, especially for low-redshift (z < 0.3). Although we 
have noted the diminishing ability to distinguish energy 
condition fulfillment from violation at the high-z end, we 
nevertheless find that our result for SEC hints at a re- 
cent transition from deceleration (Fsec{ z ) = l{ z ) > 0) 
to acceleration (Fsec( z ) = q( z ) < 0) with the transition 
rcdshift z m 0.5, if we ignore the probably biased result at 
the high rcdshift end. This estimation of transition red- 
sh ift is consistent w ith t he transition redsh i fts rep orted 
bv lRiess et ail (|2004D and lShapiro fc Turner! (|2006li using 
a kinematic model with q(z) linear in z. 

For all energy conditions considered, our result s show 

qualita ti ye agreement wi th pre vious studies rtVisserl 

19973b!; IVisser fc Barcelol l2000t ISchuecker et al.1 120031: 
Cattoen fe Visserl I2008D . In particular, the clear trend 
of DEC f ulfillment at low red shift is similar to what was 
found bv lLima et all <|2008alTbh . 

These results are obtained without reference to any 
specific form of dark energy. With this in mind, we can 
use these results to assess the feasibility of theoretical 
dark energy models. 

Some models predict that the dark energy compo- 
nent undergoes a transition into the "phantom" phase 
with the equation of state parameter w < — 1 (jStefancid 
120051 and references therein). It has been noted that 
NEC (hence necessarily DEC) violation is a feature of 
these models, and some oscillating universe models prc- 
dict episodes of NEC violation throughout the history 
(jCai et al.ll2010( ) . Even if cert ain phenomenolo gical mod- 
els could offer a fairly good fit (jXia et alJl2005l ). our data- 
driven result for NEC suggests that there is not enough 
motivation to consider them as a viable explanation for 
cosmic acceleration, at least not in the redshift range 
investigated in this paper. 

On the other hand, a family of dark energy mod- 
els exist that do not violate NEC or DEC, yet display 
transient periods of SEC violation (iBarrow et al.l 120001 : 
iBlais fc Polarskil 120041: iBarenboim et al.l 120061 ) While 
such scenarios are unlikely to have been realized from the 
early universe up to z > 1 (|Linder 201(1 iLinder fc Smith! 
l20Tl . from an energy-condition point of view they are 



not ruled out by our independent results in low rcdshift. 
Whether the "returning" -Fsec curve in Figure [7] within 
z < 0.4 suggests a coming back to SEC fu lfillment, which 
has been independently noted earlier (jShafieloo et al.1 
120091 : iGuimaraes fc Limall201l! ) for spatially fiat models, 
may be worth examining using future data. 

As we have already noted in Section 0] our test is sen- 
sitive in low redshift but not in high redshift. In the 
future, this could be partially alleviated by the availabil- 
ity of high-quality cosmological data that could populate 
higher redshift. 

Nevertheless, our work addresses several deficiencies in 
previous literature and has some special merits. 

First, by proposing a set of new models based on our in- 
dication functions, we are able to visualize the evolution 
of the Universe as a history of energy condition violation 
or fulfillment in a straightforward manner. By construc- 
tion, we can avoid ambiguous statements such as "SEC 
implies H (zf 1 5) < 58 ±7 k m s" 1 Mpc" 1 " (para- 
phrasing Eq. [79] of lVisserlll997bO . where the inequality 
nature of the energy conditions find some difficulty being 
used with error bounds from o bservational data. Com- 
pared with the research done bv lLima et al.l (|2008a| ). our 
models and their presentations can be more readily inter- 
preted as answers of "whether an energy condition was 
violated for a given epoch" , because our indication func- 
tions stay constant in each rcdshift bin. 

Second, our models are general and are free from 
overly restrictive presumptions on the nature of the 
da rk energy. It st a nds i n contrast to previous studies 
by ISchuecker et ail (|2003f ) where a constant equation of 
state parameter model was used to study the possibility 
of energy condition violation. 

Third, by incorporating the Fisher matrix formalism 
in our analysis, we find an optimal set of bases for the 
analysis of energy condition violation in the sense that 
the coefficients for different modes are nearly uncorre- 
cted. When we generate this set of bases, the effects 
of Ok are naturally and explicitly t aken into account. 
Comp ared with the parametrization of lShapiro fc Turnerl 
(l200ffl . our method is more generally applicable because 
we do not presume spatial flatness, and has a unified de- 
scription for all the energy conditions considered. We 
apply this formalism to the combination of SNIa and 
H(z) data, and this can be extended to other suitable 
cosmological probes when the data become available. 
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